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Abstract. For Petrov D vacuum spaces, two simple identities are rederived and some 
new identities are obtained, in a manageable form, by a systematic and transparent 
analysis using the GHP formalism. This gives a complete involutive set of tables for 
^ ■ the four GHP derivatives on each of the four GHP spin cocfhcicnts and the one Weyl 

, tensor component. It follows directly from these results that the theoretical upper 

bound on the order of covariant differentiation of the Riemann tensor required for a 
Karlhede classification of these spaces is reduced to two. 
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' 1. Introduction 
H , 

, In his pioneering work in the NP formalism, which provided definitive results on vacuum 



Type D spacetimes [18], surprisingly Kinnersley did not exploit, and in fact seems to 
have been unaware of, two simple NP identities for these spacetimes 

Ii = nn ~ TT = = pp, — pfi = I2 (1) 

(In the usual NP notation ^21j, the two tetrad vectors l,n are principal null vectors of 
the Weyl tensor.) More surprisingly, it seems that it took longer than another decade 
before these identities were noticed; and, even more surprisingly, they were obtained by 
chance — for the vacuum case. Debever and McLenaghan pj had discovered, as the 
result of a simple calculation, for Type D electrovac solutions that 

$(7r7f - rf ) = = ^{pp, - pp) (2) 

X Present address: Ghent University, Department of Mathematical Analysis. Galglaan 2, 9000 Ghent, 
Belgium. 
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from which ([T]) immediately follows — but only for $ 7^ 0; this led them to investigate 
if these identities were also true for the vacuum case, and they confirmed the validity of 
identities ([I]) for the vacuum case — by a direct calculation using the explicit metrics 
given by Kinnersley. Subsequently, Czapor and McLenaghan ^ established ([1]) for 
Type D vacuum spacetimes — without integrating the equations, but by a much more 
complicated and much longer calculation than for the electrovac case; the calculation 
was only achieved with computer support. Later these authors used this calculation 
as an example of their Maple Package for calculations in the NP formalism [7]. They 
pointed out that the presentation of some of the intermediate expressions used in their 
calculations required pages, and that some of the steps 'would be virtually impossible 
by hand'; furthermore, they stressed that 'even aided by computer, this calculation 
is far from automatic'. Their calculation was built on the unusual tetrad gauge choice 
TT = p, using a programme involving repeated applications of the commutator equations, 
judicial substitutions and tricky factorisations. 

However, the identities ([T]) are clearly invariant under spin and boost gauge choices, and 
so are valid in any gauge; the GHP formalism [H], [15], [9] which is invariant under 
such gauge choices, should be a more efficient tool, and it is to be expected that the 
results should come out without the need for such long calculations, and indeed more 
transparently than in the NP formalism. The GHP formalism has proved very efficient 
in a number of calculations involving Type D spaces, [IE], [22], [TD], [II],[2]; indeed, in 
the latter, Carminati and Vu have re-established these identities ([1]) as an example of 
their Maple Package for calculations in the GHP formalism. 

These identities do not seem to be well known; for example, an investigation of the 
invariant classification of vacuum Type D metrics in [H |5] did not use them in a situation 
where they would have been particularly useful. So the first purpose of this paper is 
to highlight these identities ([1]), as well as to draw attention to some other identities 
for spin coefficients, some of which we believe to be new. The second purpose is to 
emphasise the efficiency of the GHP formalism by obtaining all of these identities in a 
comparatively brief, straightforward, and systematic manner within the GHP formalism 
— by hand for transparency, but with computer support for accuracy. The third purpose 
is to reconcile the theoretical and measured upper bounds in the invariant classification 
for Type D vacuum spaces — without explicitly integrating the equations. 
A procedure, for the invariant classification of spacetimes, which originated with Cartan 
[3], was developed by Brans [Ij, and refined by Karlhede [E], has established theoretical 
upper bounds on the order of the Cartan scalars — the Riemann tensor and its covariant 
derivatives, calculated in a particular frame. At first sight, for some Petrov types, these 
upper bounds can be as high as seven, so it is of great practical significance to find 
whether these values are sharp, or whether they can be lowered to a more manageable 
number. In [H [5j it was shown, for Type D vacuum spacetimes, how the GHP operators 
could replace the covariant derivatives in the investigation of the different orders of 
derivatives of the Weyl tensor. By this means, the GHP formalism was exploited to 
lower the theoretical upper bound to three in the invariant classification of Type D 
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spacetimes — without explicitly integrating the equations; more precisely, it was shown 
in [H [5] that the theoretical upper bound was two for two subcases, while for a third 
subcase the bound was shown to be three, but a fourth subcase of Petrov D vacuum 
spacetimes was overlookec|§]. Petrov type D vacuum is the only Petrov subclass where all 
spacetimes are known explicitly, and the upper bound has been calculated directly from 
these metrics to be two [21], [25]. In this paper we shall strengthen the result in [51|5], 
and reduce the upper bound to two directly for the remaining subcases, by a further 
exploitation of the GHP formalism, but without explicitly integrating the equations; 
so the theoretical and measured upper bounds are now in agreement. To achieve this, 
we exploit a series of identities for vacuum Petrov D spaces, culminating in a complete 
set of tables for each spin coefficient and the one Weyl tensor component, "^2- These 
complete tables equate the four GHP derivative operators of each of the spin coefficients 
and \l/2 to expressions built only from the spin coefficients and "^2- 
The paper is organised as follows: the complete tables for the spin coefficients and 
\E'2 are computed in section O In section [3] we build on the work of [H |5] and show 
how these results enable us to derive that the Karlhede bound for the vacuum Petrov 
type D spacetimes is two, thus improving the result in [H [5]. Our conclusions and a 
brief discussion are given in the final section. The appendix contains a classification of 
vacuum Petrov type D spacetimes adapted to our purposes. 

The computer calculations alluded to above have been performed with xAct [I9]. This 
is a suite of Mathematica packages with many features, including the canonicalisation of 
tensor expressions by means of powerful algorithms based on permutation group theory 
|20] . and the ability of working simultaneously with several different frames. 

2. Identities and complete tables 

We follow the conventions and notation of [H] in our presentation. We will first write 
down the basic equations — the Ricci, Bianchi and commutator equations — of the GHP 
formalism particularised, in the usual way, for vacuum, type D spacetimes. Because of 
the simple structure of these spacetimes in the GHP formalism, all equations have ', * 
and *' counterparts (some of which may not be independent) but we shall not use * and 
*' explicitly in the presentation, although they are very useful for checking. 

Nonvanishing weighted spin coefficients and curvature components. 



5 = {p, 



Ricci equations. 



dp = t{p - p) 



(3) 
(4) 



§ The investigation in [4, 5J relied on a result in [12] which claimed that there were only two additional 
subcases in addition to the generic subcase pp'rr' ^ 0; this ignored the null orbit possibility p 7^ 0, p' = 
[6], which yields a subset of the Kinnersley Case HE metrics, which we call Class II in our Appendix. 
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Pr = p(r-f'), (5) 

5r = r^, (6) 

P'p = pp' - rf - ^2 + d'r, (7) 

and their ' counterparts; 

Bianchi equations. 

p^2 = 3p*2, (8) 

= 3r^2, (9) 
and their ' counterparts; 
Commutator equations. 

[PV - P'P]Vpg = ((r - r')d + (r - f')d' - p(^2 - rr') - q{^2 - ff')')m) 

[Pd - dP]r]p, = (^pd - f'P + gpf ') r]p„ (11) 

[Pc?' - d'P]rjp, = (^pd' - r'P + ppr')r/p„ (12) 

[P'c? - dP']n,, = (^p'd - rP' - pp'r)vp„ (13) 
[P'd' - d'P']vp, = (^p'd' - fP' - gp'f ) (14) 

[dd' - d'd]vp, = ((p' - p')P + (p - p)P' + Pi^2 + pp) - qi^2 + pp'))m) 

where we have written out all commutators explicitly acting on a scalar rjpq of spin and 
boost weight {p — q)/2, {j> + q)/2 respectively. 

When we substitute the Bianchi equations into the commutators for \&2 and also make 
use of the Ricci equations, in some cases we obtain the trivial identity; from the other 
commutators we obtain the 

Post-Bianchi equations. 

^2{-P'p + Pp -Tf + t't) = 0, 
^2{d'p - Pr') = 0, 

and their ' counterparts; and since \E'2 7^ 0, it follows that 

p'p = pp' -Tf + r'f', (16) 

d'p = Pr', (17) 

and their ' counterparts. Both these equations are quoted in p], but the first (fT6l) is 
overlooked in [H [5] . 

Applying the commutators to the spin coefficients given in ([3]), (j4]), ([5]), ([6]), yields no 
new information. 

Note that although we now have equations with terms involving all four operators acting 
on p, r (and hence also, by ' symmetry, on p', r') these are not all given in separate 
equations. Therefore it will be convenient to introduce new symbols; we follow the 
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notation of Carminati and Vu [2] and introduce the zero-weighted Zx and the {0, 2}- 
weighted Z2, 

Zx = d'r, 

= V'p - pp +Tf + ^2, 
= Pp' - pp' + t't' + ^2, 

= dr' + pp'-pp', (18) 

and 

Z2 = d'p, 

= Vt', (19) 

where we have used the Ricci equations and post-Bianchi equations to obtain the 
different alternatives. 

Note that, because of the Ricci and post-Bianchi equations, 

Z', = Z,+p'p-p'p, (20) 

whereas 

^2 = dp', 

= Vt. (21) 

We will use —Z2 in place of the Z^ used in |[2j. Using the Sachs symmetry * we have the 
properties 

Z2*' = Z2, (22) 

Zi* = Zi + pp'-rr-^2. (23) 

We now have a set of equations involving all four operators acting on p, r and by ' 
symmetry, on p', r'. These are the Ricci equations, the post-Bianchi equations and 
( |T8|) -(fT9l). Thus we can now apply each of the commutator equations to each of the 
four spin coefficients. In a number of cases the commutators are trivially satisfied via 
the Ricci equations, Bianchi equations and post-Bianchi equations. The remaining non- 
trivial results from the commutators are a set of first-order equations for the new symbols 
Zi, Z2: 

VZi = p{2Zi -Zi + t't' - pp' + p'p) + (r - f')Z2, (24) 
dZi = r(2Zi + ^i2 + ^>2 + pp + pp) - (p - p)^2 - r'{pp' - pp'), (25) 
and their ' counterparts; together with 
VZ2 = SpZ2, 

dZ2 = TZ2 + 2(p - p)Zi + p^2 - P^2 + 2pp'(p - p), 
d'Z2 = 3t'Z2, 

V'Z2 = p'Z2 + ^2r' + *2(r - 2r') - 2(f - r')(^i + pp'), (26) 
and their ' counterparts. 
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At this stage we have a complete and involutive set of tables for the eight complex 
quantities in S~^^ = {p, r, p', r', "^2, Zi, Zo, Z'2}. This enables us to rederive and complete 
the result in [H [5]. All the elementdj]] of S~^~^ are present at the second order of 
differentiation of the Weyl tensor, and there are no new quantities introduced at the 
third order; moreover, there is also no change in the invariance group (including the 
case overlooked in [U [5], which corresponds to Class II in the Appendix) at the third 
order, and therefore the Karlhede algorithm terminates at third order. 

However, in the following, by continuing the systematic analysis, we shall obtain more 
identities, and then be able to strengthen this result in the next section. 

The next step is to check the integrabihty conditions for (!2^ - (l25l) : i.e., apply each of the 
commutator equations to Zi and substitute the appropriate expressions from (!24l) - (!25l) . 
From the non-trivial results we obtain the new equation 

^2irf~T'f' + pp' -pp') = 0, (27) 

from which, by taking real and imaginary parts on the second factor, we get the (GHP 
version of the) identities ([1]): 

Ji = r'f' -Tf = = pp' - pp' = I2. (28) 

Note that these identities are invariant under the ' symmetry (and also the * symmetry). 
These identities cause simplification in some of the post-Bianchi equations (|T6l) . and in 
some of the definitions f|T8|) and fl20l) (where Zi = Z[)] we will make frequent use of 
these simplifying identities in the following. 

In the same way we check the integrabihty conditions for (!26|) : from the non-trivial 
results we obtain the new equation 

^2(^2 + pr - 2pr' + r'p) =0, (29) 
from which we obtain Z2, 

Z2 = -pf + 2pr' - t'p (30) 
and hence from (fT9l) the four new identities: 

d'p = —pf + 2pr' — t'p, 

Pt' = -pf + 2pr' - t'p, (31) 
together with their ' counterparts. 

The results so far obtained duplicate those of Carminati and Vu p]; but we shall now 
continue the analysis even further. 

For consistency it is now necessary to substitute the value for Z2 into (126|) : making use 
of the Ricci equations, the post-Bianchi equations and the identities (j28l) . yields the 

II As noted above, in [A, 5j, the first of the post-Bianchi equations (|16p was overlooked; this is why 
in that paper the corresponding expressions for D^'^ and D^"^ contain nine variables: there is an 
additional variable (Pp'), which corresponds to Z'l in our notation, but which we know is related to Zi 
because of (fT6| . 
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(32) 

(33) 



together with the ' counterpart of ( l32l) (which is simply its complex conjugate), and the 
' counterpart of (!33|) given by 



Note that (!33|) and (jUj) are not independent when p' 7^ 7^ p, but in the cases p' = 
or p = they may provide different information. 

In addition, it is also necessary to apply each of the four operators to the identities fl28l) 
(in practice, only one of the operators needs to be checked because of symmetry); the 
result is a duplication of the existing information. 

At this stage we have a complete and involutive set of tables for the six complex 
quantities in S'^ = {p,T, p' ,t' ,"^2, Zi}. However, this does not permit us yet to lower 
the Karlhede bound; in order to do that we need to be able to substitute Zi in terms 
exclusively of the elements in S. 

Until now, we have been looking at the whole class of Petrov D vacuum spaces together: 
for we have a complete table of its derivatives, i.e., expressions for all four of its GHP 
derivatives; but for each spin coefficient we have incomplete tables (from the point of 
view of elements of S), since we have obtained explicit expressions built from S, for 
only three of its derivatives. These expressions, contained in ([3])-([6]), fl3T|) . involve only 
algebraic products of spin coefficients, and for the sake of clarity are gathered in the 
next expression 



From here we also get the derivatives of p', r' by taking the prime of each equation. 
There remains one derivative for each spin coefficient which contains terms involving 
Zi, in the generic case (see Class IIIB, below), we can obtain Zi from ( I32l) -( l34|) . and it 
is found to consist of only algebraic products of spin coefficients and ^^2, but there are 
cases in which (|32l) - (!Ml) are not independent, and hence we must proceed in a different 
way. To continue the analysis requires that we now separate into the different classes 
which are presented in the Appendix. 

Class I: pp' ^ 0; r = = r'. 

Equation (fTSl) implies Zi = and hence adding (fTSl) to (j35l) it follows immediately 
that all the GHP derivatives of all non-zero elements of S are given by expressions 
constructed algebraically from elements of S: a complete involutive set of tables for the 
action of the GHP operators on all of the nonvanishing elements of S. 
Note that in this case we also have an extra relation involving the spin coefficients which 
is obtained by setting Zi = in ( |33l) . 



pZi - pZi = 2pp(p' - p') + rr(p - p) + p^2 - P*2, 



(34) 



Pp = p'^, dp = T{p~p), d'p = -pT + 2pT' -t'p, 
VT = p{T-f'), V't = -p'r + 2p'T - rp', dr = r 



.2. 



(35) 



^2P' + 2(p' - p')pp' - ^2P' = 0, 



(36) 
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together with its ' counterpart which is obtained from fl34l) . However, the latter is not 
independent of fl36|) . since we assume p' 7^ 7^ p in this case. 

The integrability conditions arising from ( l35l) give a set of conditions which depends 
functionally on (136|1 . Similarly if we differentiate (136!) and replace the derivatives of the 
spin coefficients and by their expressions we get a condition which is functionally 
dependent on fl36l) and thus gives nothing new. 

Class 11: tt' ^0; p' = ^ p. 

The third hue of 1^ gives Zi = ^2 + r'f' and hence adding dUl) to ([35]) it follows 
immediately that all the GHP derivatives of all non-zero elements of S are given by 
expressions constructed algebraically from elements of iS: a complete involutive set of 
tables for the action of the GHP operators on all of the nonvanishing elements of S. 
Also (l32l) becomes 



^2T + 2(r - r )rr - = 0, (37) 

but fl33|) and flM|) are identically satisfied. The substitution of the values of Zi and Z2 
in (!24l)-( !25|) yields the extra restriction (since p 7^ 0) 

p(^2 + ^2 - 2rr' + 2rf - 2f'f ) = 0. (38) 

The last pair of equations can be regarded as a linear system in the variables "$2, ^2 
from which it can be deduced that 



(r + r)^2 = 2rV (39) 



and so 



2rV 

^2 = -— , (40) 
r + r' 

where, from ([39]), r' 7^ 0, r 7^ imply that r + f' 7^ 0. 

Next we insert (llOl) into ([8])- ([9]) and their / counterparts, and use the expressions already 
found for the derivatives of the spin coefficients. The result is the restriction 

h=pf + pr' = 0. (41) 

Neither the integrability conditions of (135]) nor the differentiation of (14T!) give further 
conditions. Finally we note that the combination of (HOl) and (HTi) yields the relation 

^2P' + ^'2P' = 0. (42) 
Class II ': rr' ^ 0; p = ^ p'. 

This class does not need to be considered separately, since, by the interchange of the 
null vectors / ^ n, this class is transformed into the previous one. 

Class IIIA: tt' p' p 7^ 0; pV - p'f ' = = pr' - pf. 

When we differentiate the conditions p'r — p'f' = = pr' — pf we obtain 

p = p, p' = p', r = f', ^2 = ^2, ^1 = ^1, (43) 
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We take as independent variables p, p', r, and \E'2. Combining with (1291) and ([2SD 
gives respectively 

^2 = = Z^, (44) 

and 

Pr = = PV, dp = = dp', d'p' = = d'p. (45) 

The only remaining equations where we can get information about Zi are (1241) - (l25l) . 
which simplify to 

m = p{Z, + Tf), (46) 
dZi = 2T{Zi + ^2 + PP')- (47) 

Putting 

3A = (Zi + Tf)/^l^\ 3B = ^2^' + {Zi + pp')/^2^', (48) 

we find directly from the above equations that 

PA = = P'A, VB = = P'B, (49) 

dA = 2rfi^2^V3, d'A = 2t'B'^1^^/3, dB = 3r^2^^ d'B = Sr'^l^^. (50) 
From the last set of equations it follows that 

d{B^-9A) = = d'iB''-9A). (51) 
For zero- weighted quantities, such as A,B, we have the relationship 

Va = riaP + laP' " rhaO - rUad' , 

and so this implies, 

Va(52 - 9A) = {B^ - 9A) = k, real constant. (52) 

Therefore 

(Zi + ^2 + PpT - HZ, + rf )^2 - k^T = 0, (53) 

and so Zi can be obtained explicitly in terms of the spin coefficients and ^^2 (and the 
constant k). 

Therefore, by substitution for Zi from (!53l) into (ITSl) . and then adding (ITSl) to (135!) . it 
follows immediately that a// i/ie Gi/P derivatives of all non-zero elements of S are given 
by expressions constructed from elements of S: a complete involutive set of tables for 
the action of the GHP operators on all of the nonvanishing elements of S. 
The integrability conditions of this set of equations are identically fulfilled and thus they 
give no further restrictions. 

Classes IIIB,C: tt' p' p ^ 0; p'r - p'f' ^ Q ^ pr' - pf. 
In this case from (!32l) - (l33l) it follows that 

(p'f' - p't)Z, = t(^^2P' + 2(p' - p')p'p - ^2P' + (2p'r - p'r - p'f')f):54) 
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Given that p'r — p'f' ^ this equation enables us to obtain Zi, and by inserting such 
a result in f|T8l) . we obtain the expressions for the one remaining derivative of each spin 
coefficient. These expressions are 



together with their ' counterparts. 

Hence, adding fl55|) to fl35|) . it follows immediately that all the GHP derivatives of all 
non-zero elements of S are given by expressions constructed algebraically from elements 
of S: a complete involutive set of tables for the action of the GHP operators on all of 
the nonvanishing elements of S. 

The integrability conditions of this set of equations are identically fulfilled and thus they 
give no further restrictions. 

Class IV: rr' ^ 0; p' = = p. 

The results from Class 11 are still valid and we only need to set p = in them; therefore 
all the GHP derivatives of all non-zero elements ofS are given by expressions constructed 
algebraically from elements of S: a complete involutive set of tables for the action of the 
GHP operators on the non-vanishing spin coefficients. 

In this case, ( l38l) does not arise, and hence we only need to consider f l37|) . Again, neither 
the differentiation of this condition nor the integrability conditions of fl35l) result in new 
relations. 

3. Invariant classification 

The Karlhede algorithm [1^ for an invariant classification of geometries involves 
calculating — in a particular frame — the (zeroth derivative of) Riemann tensor and 
then computing successively higher orders of its frame covariant derivatives ('Cartan 
scalars'), together with possible changes of frame. The algorithm terminates at that 
particular order where (i) no new functionally independent Cartan scalar is supplied by 
the frame derivatives of the Riemann tensor; and (ii) no further fixing of the frame 
is possible using these frame derivatives, i.e., the invariance group of the frame is 
unchanged. This particular order is called the Karlhede bound. 

Collins et al [U [5] have used the GHP formalism to establish a relationship between 
the Cartan scalars and the GHP derivatives of the Weyl tensor, and hence provided an 
efficient approach to the Karlhede classification of vacuum Petrov type D spacetimes. 
We shall now exploit their results directly, and strengthen their final conclusion. It is 
shown in [ll[5], and in their notation, that: 

• the zeroth-order covariant derivative of the Weyl tensor is given by Z^°\1'(\1'2); 




d'r 




(55) 



p't' — p't 



Petrov D vacuum spaces revisited: Identities and Invariant Classification 11 

• the first covariant derivatives of \E'2 are shown to be — via the Bianchi equations 
— algebraic expressions involving only all elements of S, and so are given by 
D^'^{'^2, P, p' ,t,t'); for each of the different classes, the dimension of the invariance 
subgroup changes. 

This means that we have to continue to the next order of differentiation. This is where 
we strengthen the earlier result because our analysis in the last section shows that 
derivatives of the spin coefficients do not in fact introduce any new quantities (since we 
have found explicit expressions for Zi exclusively in terms of the elements of 5, in all 
cases.) Hence, 

• the second covariant derivatives of — via Bianchi equations, Ricci equations, fl35l) 
and ( fT8|) substituted with the various explicit expressions for Zi — are expressions 
involving exclusively elements of S and are given by Z)^\E'(\1'2, p, p', t, t'); also the 
dimension of the invariance subgroup does not change for any of the different classes. 
(For completeness, we note that in it was pointed out explicitly — for the cases 
corresponding to our Classes I and IV — that since all components of D^^^ can be 
expressed algebraically in terms of first order quantities, the frame cannot be fixed any 
further at second order; in addition — for the case corresponding to our Class III — the 
frame was fixed completely at first order, and so no further fixing is possible at second 
order. For our Class II — which was overlooked in [H [5] — the frame can be fixed 
completely at first order by choosing, e.g., |p| = 1 and r = f; and so no further fixing is 
possible.) 

In summary, all elements of S occur at first derivative order, and no new quantities 
appear at second derivative order; moreover, the invariance group remains unchanged 
from first to second derivative order: therefore the Karlhede algorithm terminates at 
second order, 

4. Conclusion and Discussion 

Thanks to the efficiency of the GHP formalism, we have rederived the identities (|28|) . 
and derived some new identities ( !32l) . ( !33|) . ( !34l) . ( l55l) for the spin coefficients, in a 
transparent manner; using these results it has been easy to lower to two the upper 
bound for the invariant classification of vacuum Petrov D spaces. 

In addition to showing in Section 2 that all GHP derivatives of all non-zero elements 
in S are given by expressions built exclusively from elements of S, we have obtained a 
complete and involutive set of explicit tables for the action of all four GHP derivatives 
on the elements of iS — in a manageable form, by a systematic and transparent analysis 
using the GHP formalism. These identities (!28l) . as well as partial tables (!35l) (for three 
GHP derivative operators acting on each spin coefficient) have already been obtained 
by Carminati and Vu [2] in the GHP formalism, but we believe that this is the first time 
that the spin coefficient identities ( l55l) are presented and that all possible compatibility 
and integrability conditions have been explicitly checked for these spaces. These tables 
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will be the basis for a systematic and efficient integration of the different type D vacuum 
classes in a future paper. 

In their NP Maple investigation [7], Czapor and McLenaghan state that they have 
obtained the complete tables for all of the NP spin coefficients (including the 'badly- 
behaved' ones), but that these tables, together with a lot of the intermediate 
calculations, were too long to print out explicitly in the NP formalism. Therefore, in 
their work [7j, there was perhaps the potential to deduce the result that the theoretical 
upper bound was two in the invariant classification (at least, in the generic case); 
however, this possibility was obscured by the presence of the badly behaved spin 
coefficients which led to the very long and unmanageable expressions, and checking for 
possible tetrad changes would probably have been virtually impossible. In comparison, 
the GHP calculations and results are manageable and transparent, and it is very 
instructive to compare how the calculations develop in the GHP procedure of this paper 
compared with the NP version in [7]. 

The classification given in the Appendix improves on, and ties up a few loose ends 
regarding earlier classifications. 

The system xAct was used to cross-check all the computations derived by hand and 
to derive others which due to their complication were difficult to obtain without the 
aid of the computer. A Mathematica notebook in which all the equations of the GHP 
formalism are obtained from scratch is available in [T3] . 
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Appendix: Classification of vacuum Petrov type D spacetimes 

We shall give a classification for vacuum Petrov type D spacetimes which is well suited 
to our investigations; since this classification eventually turns out to be closely related 
to the Kinnersley classification we shall label it similarly to that classification. 
In addition to the GHP zero-weighted quantities Ji and I2 we shall use the GHP versions 



h = pT' + pr, h = p'r + p't' 
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h = pr - pT, Ie = pr - pr 

of the NP quantities /s, 14,75 which have also been defined in [6], and in addition, for 
completeness, we have added Iq. 

These are not invariants under tetrad transformations, but since they are of GHP weight, 
they are rather rescaled; so, if for example J3 = 0, then under arbitrary spin and boost 
tetrad transformations, I'^ = 0. By making direct use of the results Ji = = I2, the 
simple lemma follows. 

Lemma 1. If p ^ 0, then: 

(i) h = h = 0; 

(ii) h = h = 0. 
If p' ^0, then: 

(i) h = h = 0; 

(ii) h = ^ h = 0- 

From the earlier work, we obtain two other simple lemmas. 
Lemma 2. 

(i) p = ^t^O^t', (56) 

(ii) p' = ^ r ^ ^ r'. (57) 



Proof: Assume that p = 0; then (JTj) becomes 

= -Tf - ^2 + O't, 

which implies r 7^ in a type D spacetime. Combining this with the condition /i = 
shown in ( l28l) we conclude that r' 7^ 0. If p' = the argument is similar. □ 

Lemma 3. 

(i) p ^ ^ p', r = r' = 

(ii) p^0^p',r' = T = 

Proof: Apply (jS]) and its prime counterpart. □ 

From these lemmas it follows that any vacuum type D spacetime falls within one of the 
following mutually exclusive subclasses: 

Class I. pp' ^ 0; r = = r'. 

Class Ilfl rr' ^ 0; p ^ = p'. 

[ Class II'. tt' ^ 0; p' ^ = p. ] 

Class III. pp'tt' ^ 0. 

Class IV. rr' 7^ 0; p = = p'. 

^ This was the case overlooked in [12] and [J, [5]. 
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Class 111 can be further subdivided: 
Class IIIA. Jg = = J5; /4 ^ ^ I3. 
Class IIIB. Iq^O^I^; h^O^ J3. 
Class inc. Je ^ ^ Is; h = = h- 

We have used these labels because of the substantial overlap with the Kinnersley cases: 
by examining the different explicit metrics obtained by Kinnersley, we deduce the 
following: 

• Class 1 coincides with Kinnersley Case 1. 

• Class 11 coincides with a subcase (p' = 0) of Kinnersley Case HE. 

• Class lllA coincides with Kinnersley Case lllA. 

• Class lllB coincides with Kinnersley Case lllB. 

• Class lllC coincides with Kinnersley llA-D, F, together with the other subcase 
(p' 7^ 0) of Kinnersley Case HE. 

• Class IV coincides with Kinnersley Case IV. 

A partial classification of the vacuum Petrov D metrics with non-zero cosmological 
constant, using Ji, I2 ■ ■ ■ h, was given in [6], and the above results are more complete. 
Note that the table on p2166 of [6j permits p' = 7^ I3, which is ruled out above; and 
so the only Kinnersley case with null orbits (in the vacuum case with zero cosmological 
constant) is HE. our attention has been drawn to a classification given in [23] which is 
very similar to the one in this paper, and it also notes these additional properties. 
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